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ARTICLE INFO ABSTRACT 


Keywords: The quality and safety of Pu-erh tea are largely determined by microorganisms, processing and storage condi- 
Pu-erh _ tions. The systemic co-relationship between microbiome and metabolomics during manufacturing of Pu-erh tea is 
Microbiome not fully clarified. We comprehensively explored the microbial and metabolic dynamic evolution and revealed 
rite the decisive factors for the quality and safety of Pu-erh tea based on creative methods. The results showed that 


fungal abundance increased in the middle of pile fermentation period and then decreased during storage, but the 
bacterial abundance increased in the initial stage, subsequently decreased and finally increased. All microor- 
ganism species as well as dominant fungi and bacteria were identified during different processing steps. A total of 
209 metabolites were identified and 50 of them were characterized as critical metabolites responsible for 
metabolic changes. Multivariate analysis verified that these critical metabolites were generated by specific 
dominant species. The environmental condition exhibited a great effect on the growth and reproduction of 
bacteria, but less than 10% of the dominant fungi exhibited significant correlation with temperature and hu- 
midity. In conclusion, our results revealed the microbial composition is the critical factor in changing the 


metabolic profile and provided a theoretical basis for improving the quality and safety of Pu-erh tea. 


1. Introduction 


Pu-erh tea (PET) is a typical dark tea manufactured by microbial 
fermentation of Camellia assamica tea leaves (Zhang, Zhang, Zhou, Ling, 
& Wan, 2013). Its popularity is due to the attractive flavour, taste and 
putative health benefits, including anti-obesity (Lu et al., 2019), 
anti-oxidative (Padurariu et al., 2010), anti-cancerogenic (Wu et al., 
2007) and free radical scavenging (Qian, Guan, Yang, & Li, 2008). The 
complex biochemical profile of Pu-erh tea results from microorganism 
present during the manufacturing process and is considered to be the 
critical element in forming its bioactivity, quality and sensory 
characteristics. 

Pu-erh tea is usually produced by a natural solid-state fermentation 
(SSF) process using sun dried green tea leaves as raw material. The tea 
leaves are moistened with water and matured by pile-fermentation for a 
few weeks (Lv, Zhang, Lin, & Liang, 2013). Generally, the tea is turned 
over once a week during pile fermentation and fermentation is stopped 
when the tea mass is reddish-brown and free from the astringent taste 
(Lee & Foo, 2013). This natural SSF process is associated with a series of 
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reactions, such as degradation, oxidation, condensation, structural 
modification, methylation, and glycosylation (Zhu et al., 2020). Previ- 
ous studies have reported that the content of catechin derivatives, fla- 
vonoids and their glycosides, phenolic acids, alkaloids and terpenoids 
are significantly changed during the SSF process (Zhu et al., 2020). 
However, in these studies, only a small number of common compounds 
have been characterized and few systematic studies have focused on 
determining the variation in the metabolic patterns due to limitation in 
available analytical methods. Therefore, a more comprehensive inves- 
tigation is needed to understand and clarify the dynamic evolution of 
microbiome and metabolites during the entire manufacturing process. 
To our knowledge, fermentation influences biological and chemical 
processes as well as microbial community composition. Moreover, the 
microbiome has been greatly advanced with the development of mo- 
lecular tools. Diverse microorganisms dominate in different steps of 
fermentation process affect the Pu-erh tea quality. Nevertheless, the 
dynamic evolution of microbial communities and metabolites changes 
during the Pu-erh tea fermentation process is still not completely deci- 
phered. Unravelling the correlation between microbiome and 
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metabolites is critical for developing safe, consistent and high-quality 
Pu-erh tea. High throughput sequencing is a powerful tool to investi- 
gate the abundance and identity of microorganisms (Robin, Ludlow, 
LaRanger, Wright, & Shay, 2016). Furthermore, quantitative PCR 
(qPCR) is an accurate method for the detection and quantification of 
microorganisms due to its high sensitivity, specificity, rapidity and di- 
versity (Postollec, Falentin, Pavan, Combrisson, & Sohier, 2011). To 
accurately elucidate the comprehensive dynamic evolution of microor- 
ganism during the processing and storage of Pu-erh tea, we combined 
high throughput sequencing and qPCR to directly and precisely quantify 
copy numbers of the microbial genome for different samples. In our 
previous study, we revealed that microbial composition is the critical 
factor in changing the metabolic profile of Fu brick tea using as similar 
approach (Li et al., 2021). 

The metabolomics strategy can unbiasedly identify and quantify as 
many of the compounds present in a biological sample as possible 
(Want, Nordstrom, Morita, & Siuzdak, 2007) and allow us gain a better 
insight into biochemical changes occurring with time and processing 
(Capanoglu, Beekwilder, Boyacioglu, Hall, & De, 2008). Metabolomics 
has been applied to measure the composition of Pu-erh tea and its 
composition changes in different years after fermentation (Ku, Kim, 
Park, Liu, & Lee, 2010) and to compare the chemical constituents of 
Pu-erh tea with black tea and green tea (Xie et al., 2009). A LC-MS 
approach has been broadly applied in metabolomics due to its high 
throughput capacity, high resolution, powerful separation ability and 
strong sensitivity (Hda et al., 2020). It mainly focuses on the identifi- 
cation of small molecular endogenous metabolites in organisms (Patti, 
Yanes, & Siuzdak, 2013). However, the nutritional and healthy in- 
gredients in Pu-erh tea are mostly macromolecules, such as poly- 
saccharide, tea polyphenols and catechin (Cheng et al., 2020). To make 
up this gap, biochemical methods were used to measure the contents of 
macromolecules in Pu-erh tea. The multivariate correlation analysis of 
metabolomics in Pu-erh tea research show that the metabolomics 
approach is a powerful tool to explore the biochemical changes occur- 
ring during pile fermentation processing. 

In this study, raw materials, five steps in the pile fermentation pro- 
cess, and two ages of ripened tea were sampled to investigate micro- 
biome composition change, correlation between microbiome and 
metabolites, correlation of environmental conditions and microbiome 
and antioxidant changes during Pu-erh tea processing. A potent 
metabolomics tool was combined with biochemical measurements to 
characterize changes in hundreds of endogenous metabolites during the 
fermentation process. We applied microbiomics (high throughput 
sequencing combined with qPCR) combined with multivariate analysis 
to elucidate the metabolic potential of the microbial community and 
discussed the function compounds changes during manufacturing pro- 
cess. Furthermore, the correlation between environment and microbial 
make-up was investigated to evaluate the complex relationship between 
functional microorganisms, metabolic pathways and dominant func- 
tional compounds in Pu-erh tea. 


2. Experimental section 
2.1. Materials and chemicals 


Eight samples from steps in Pu-erh tea processing or ripening were 
obtained from Yunnan LongRun Tea Industry Group (China). The sam- 
ples were labeled as raw materials Pu-erh tea (RMT), pile-fermentation 
at first time (PF1st), pile-fermentation at second time (PF2nd), pile- 
fermentation at third time (PF3rd), pile-fermentation at fourth time 
(PF4th), pile-fermentation at fifth time (PF5th), ripened tea 180 days 
(RT180d) and ripened tea 360 days (RT360d). 

All standard reagents involved in the experiments were purchased 
from Yuanye Biotechnology Company (Shanghai, China). The LC-MS 
grade acetonitrile, formic acid and methanol were purchased from 
Sigma Aldrich (St. Louis, MO, USA). Other reagents were analytical 
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grade from China National Medicines Corporation Ltd. 


2.2. DNA extraction, polymerase chain reaction amplification, and 
sequence 


Microbial DNA was extracted from the samples using a FastDNA® 
Spin Kit (MP Biomedicals, Norcross, GA, USA) according to the manu- 
facturer’s instructions. The V3-V4 region of the 16S rRNA gene was 
amplified from the extracted DNA using forward primer 338F (5’- 
ACTCCTACGG GAGGCAGCA-3’) (Adams, Miletto, Taylor, & Bruns, 
2013) and reverse primer 806R (5’-GGACTACHVGGG TWT-CTAAT-3’) 
(Xu, Tan, Wang, & Gai, 2016). The quality of DNA samples was evalu- 
ated by electrophoresis on 2% agarose gels. Each DNA sample was 
amplified in triplicate by PCR with a total volume of 20 tL. The PCR 
products were purified with AxyPrep DNA Gel Extraction Kit (Axygen 
Biosciences, Union City, CA, USA). Purified PCR products were quanti- 
fied by using QuantiFluotTM DNA assay Kit (Pro-mega, Madison, WI, 
USA). Following quantification, the mixed sample amplicons sequencing 
was performed using the Illumina MiSeq system by Majorbio Bio-Pharm 
Technology Co., Shanghai, China. 


2.3. Sequence processing and analysis 


The raw sequencing reads were demultiplexed using Trimmomatic 
(v1.7.0) ( Bolger, Lohse & Usadel, 2014), FLASH, (Magoc & Salzberg, 
2011) UPARSE (Dong et al., 2017) and UCHIME (Aas, Davey, & Kau- 
serud, 2017) software for denoising, trimming, quality filtering, pairing, 
and aligning. Resulting sequences clustered the operational taxonomic 
units (OTUs) at 97% similarity using UPARSE software(version 7.1 
http://drive5.com/uparse/). The representative sequence of each OUT 
was assigned to taxonomic information using the RDP classifier (http 
://rdp.cme.msu.edu/) with a 70% threshold by comparing with Silva 
database (SSU128). Community richness index (CHAO index), commu- 
nity diversity (Shannon index), and the Good’s coverage of sequencing 
were estimated using the Mothur Software. 


2.4, Standard curve generation and quantitative PCR (qPCR) detection of 
fungi and bacteria during processing and storage of Pu-erh tea 


The 188 bp fragment (091-279) of yeast (BY4741, Saccharomyces 
cerevisiae) ITS gene cloned by PCR with specific primers (F: 5’- 
CGCGGATCCCCAGCCGG GCCTGCGCT TAAG, R: 5’- 
CCGCTCGAGCCTCT GGGCCCCGATTGCTCG) was inserted into Bam HI 
and XhoI sites in the p416-TEF vector for expression in yeast. This p416- 
TEF-ITS(188) plasmid was used for fungal standard curve generation. 
The same strategy was used for bacterial standard curve generation. The 
659 bp fragment (031-690) of E.coli 16s rDNA obtained by PCR with 
specific primers (F: 5’-CGCGGATCCCGGCAGGCCTAACA CATGCAAG, 
R: 5’-CCGCTCGAGGCAT TTCACC GCTACACCTG) was inserted into 
Bam HI and Xho I sites in the p416-TEF vector. All plasmids were 
sequenced to confirm inserted gene sequences. The cycle threshold (Ct) 
value, obtained from 10-fold serial dilutions of plasmid DNA in deion- 
ized water to produce 5 different concentrations was used to generate 
standard curves for quantitative PCR (q(TOWER,3.0G) using primer sets 
that are specific for fungi (F: 5’°-CCAGCCGGGCCTGC GCTTAAG, R: 5’- 
CCTCTGGGCCCCGA TTGCTCG), bacteria (F: 5’°-ACTCCTACG GGAGG- 
CAGCAG, R: 5’-ATTACC GCGGCTGG) (Bian et al., 2021). Common 
molecular biology techniques followed previously established methods 
(Zhang et al., 2019). Genomic DNA samples of different processes of 
Pu-erh tea were extracted by DNA extraction kit (TIANamp Soil DNA kit) 
followed by quantitative PCR using same primer sets as those of stan- 
dard curve qPCR. 
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2.5. Sample preparation 


2.5.1. Determination of main active compounds 

Tea polyphenol was determined by the iron tartrate colorimetric 
method (Liang, Zhang, & Lu, 2005). The contents of theanine (TA), 
caffeine (CAF), gallic acid (GA), caffeine (CAF) and catechins, including 
(+)-catechin (C), (—)-epicatechin (EC), (—)-epigallocatechin (EGC), 
ECG, (—)-Galloca- techin gallate (GCG), (—)-Gallocatechin (GC) and 
EGCG in the Pu-erh tea samples were determined by high-performance 
liquid chromatograph (HPLC) using waters GRX-9123A HPLC System. 
Meanwhile, total flavonoid content was determined using the colori- 
metric method (SZDB/Z 349-2019) with rutin as the standard. The 
qualitative analysis of tea pigments, including theaflavins, thearubigins, 
and theabrownins, was carried out using the method of Wang et al 
(Wang, Peng, & Gong, 2011). 0.1 g tea sample was extracted with 3 mL 
methanol (70% v/v) for 20 min in 60 °C water bath. The extraction was 
repeated two times. After cooling, solution was passed through a 0.2 um 
syringe filter for HPLC analysis. The detection wavelength was 280 nm. 


2.5.2. Determination of tea free amino acid and tea protein 

Total amino acids of tea sample was quantified using ninhydrin 
method (Jabeen et al., 2019). 0.15 g tea powder was placed in a test tube 
and extracted with 10 mL hot water for 20 min in 95 °C water bath. The 
detection wavelength was 570 nm. Kjeldahl method was used for the 
quantitative determination of nitrogen contained in the tea protein. 


2.5.3. Determination of water soluble sugar 

Anthrone method was used to determine the water soluble sugar in 
tea sample (SN/T 4260-2015) with glucose as the standard. 0.2 g tea 
powder was extracted with 20 mL hot water for 60 min in 95 °C water 
bath with vibrating the tube every 10 min. The obtained supernatant 
was diluted to 50 mL and detected at 620 nm. 


2.6. UPLC-QTOF-MS analysis 


Pu-erh tea samples were determined by high-performance liquid 
chromatograph (HPLC) using waters GRX-9123A HPLC System. The 
separation was carried out using a BEH C18 column (1.7 ppm, 2.1 mm x 
100 mm, Waters, USA). The column was maintained at 45 °C eluted with 
a 1-100% methanol containing 0.1% (v/v) formic acid -aqueous con- 
taining 0.1% (v/v) formic acid gradient over 15 min at flow rate of 0.4 
mL/min. 1 pL sample was injected onto the column, and gradient elution 
condition was followed previous publication (Cheng et al., 2020). 

The mass spectrometer was operated in the negative ionization mode 
over a full scan range of 50-1000 m/z with the following settings: gas 
temperature, 450 °C; gas flow, 15 L/min; source temperature, 115 °C; 
capillary voltage, 2 kV; cone voltage, 40 V; and collision energy, 6 eV. 
The internal standard (DL-4-chlorophenylalanine) was applied to 
partially compensate for the concentration bias of a given metabolite 
among the tea samples. 


2.7. Multivariate analysis 


The raw data from LC-MS were converted to acquire format using 
MataboAnalyst software (https://www.metaboanalyst.ca/). Hierarchi- 
cal cluster analysis (HCA), unsupervised principal component analysis 
(PCA), analysis of variance (ANOVA) and supervised orthonormal par- 
tial least-squares discriminant analysis (OPLS-DA) were performed 
using SIMCA-P +software (version 14.1, Umetrics AB, Sweden). The 
HCA visualized the relationship of tea samples by simultaneously plot- 
ting the respective dendrograms and their connections. The PCA 
examined the intrinsic variation in the collected data matrix and 
distinguished the differences among these PET processing samples, with 
OPLS-DA applied to classify samples of solely Y variables. The critical 
metabolites responsible for metabolomics variation caused by PET 
processing were identified using OPLS-DA modeling combined with a 
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threshold of variable importance project (VIP)>1.16. The heat map 
analysis was carried out using the Web site https://software.broadinstit 
ute. org/morpheus/. 


2.8. Statistical analysis 


All experiments were carried out triplicates and the results were 
expressed as mean values followed by the standard deviation (n=3). 
One-way analysis of variance (ANOVA) with Dunnett’s multiple com- 
parisons test were used to distinguish the significate level of metabolites 
between groups. The linear correlation between the inhibitory effects of 
digestive enzymes and the chemical composition of the samples was 
analyzed by SMICA-P software and the probability values calculated. P- 
values below 5% were considered significant. 


3. Results and discussion 
3.1. Sequence statistics of fungi and bacteria of Pu-erh tea samples 


To comprehensively explore the microbial dynamic evolution during 
the manufacturing process and storage of Pu-erh tea, specific regions of 
the fungal ITS gene and bacterial 16s rDNA were obtained by PCR with 
specific primers and were sequenced by high throughput sequencing 
(Illumina). The total number of quantified fungal raw sequences was 
1507381 (Supplementary Table 1). By comparison, 1348266 qualified 
bacterial sequences were obtained (Supplementary Table 2). The total 
sequence lengths of fungi and bacteria were 380782600 bp and 
56434614 bp, respectively, which are more than the numbers of se- 
quences reported by previous studies (Jia et al., 2020; Ma et al., 2017; 
Ming et al., 2015). The average sequence length of bacteria was 420.1 bp 
which was higher than the average fungal sequence length of 252.6 bp. 


3.2. Comparison of fungal community and richness during the 
manufacturing process and storage 


The most abundant microorganisms in Pu-erh tea processing are 
fungi and bacteria. The quality-filtered fungal sequences of all eight Pu- 
erh tea samples were clustered into 618 fungal operational taxonomic 
units (OTUs) at the 97%-similarity cut-off. Venn diagram analysis 
showed the number of fungal OTUs was decreased during the whole 
fermentation process and ripened stage (Fig. 1A). The eight Pu-erh tea 
samples shared 3 fungal OTUs, the highest OTU number, 575, appeared 
in the RMT sample. Moreover, the Chao index suggested the RMT 
sample had the highest fungal richness compared with other Pu-erh tea 
samples (Fig. 1B). The Shannon index indicated that the raw Pu-erh tea 
sample exhibited the highest fungal diversity (Fig. 1C). The fungal 
rarefaction curves fully approach the saturation plateau (Fig. S1A). This 
result showed the fungi from RMT could not be better acclimated to 
damp and hot conditions. 

To compare the fungal communities, fungal taxonomy of different 
stages was identified. The 618 fungal OTUs could be assigned into 4 
different phyla, 18 classes, 52 orders, 115 families, 176 genera and 258 
species (Fig. 1D). Ascomycota was the predominant phylum in all sam- 
ples, accounting for 99% of total effective sequences. Basidiomycota and 
some other phyla fungi were also identified in the early processing 
stages, but Mucoromycotal appeared in the later processing stages. 

At the species level, Cladosporium and Epicoccum nigrum comprised 
>80% of all sequences in the RMT sample. Aspergillus dominated the 
whole manufacturing process of Pu-erh tea (Fig. 1E). Aspergillus niger 
could contribute to digesting the polysaccharide, lipid and protein into 
small molecules to enhance the nutritional profile or change tea flavour 
during Pu-erh tea fermentation (Abe et al., 2008; Zhang et al., 2016). In 
addition, tannase produced by Aspergillus niger hydrolyzes tannin to 
gallic acid, which contributes to the reddish-brown color of the tea 
(Aissam, Errachidi, Penninckx, Merzouki, & Benlemlih, 2005; Die- 
peningen et al., 2004). This result confirmed previous findings that 
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Fig. 1. Comparisons of fungal composition and richness in different samples of Pu-erh tea. (A)Venn diagram analysis of fungal OTU level in samples during 
manufacturing process of Pu-erh tea. (B)Fungal taxonomic compositions showing the fugal successions at species level. The fungal communities in eight period 
samples are comparably diverse. Chao index (C) and Shannon index (D). Heat map analysis of fungal (E) communities in Pu-erh tea samples at the species level. (note: 


The taxonomic abundance <1% were classified in to “others”; “*” P < 0.05 ; 


gallic acid content was much higher in the ripened Pu-erh tea samples 


(Zhang, Li, Ma, & Tu, 2011). 


Due to the consumption of macromolecules and limited nutrients, 
the abundance of Aspergillus niger significantly decreased to 53.15% in 
PF4th. On the other hand, Rasamsonia emersonii and Thermomyces 
lanuginosus increased in relative abundance to 25.78% and 4.83% in 
PF5th stage, respectively. The macromolecules hydrolyzed by Aspergillus 
niger have been proven to be a rich source for other microorganisms. 


“x? BD < 0.01; “***” P < 0.001.) 


ripened teas. Blastobotrys adeninivorans was the dominant fungal species 


in the ripened Pu-erh tea samples, followed by Rasamsonia emersonii, 


fermentation process. 


Thus, the fugal community structure shifted during the production of 


Thermomyces lanuginosus and Aspergillus penicillioides (Fig. 1D). Yeast is 
strict on the growth conditions, especially the temperature, humidity 
and nutrition sources. Therefore, yeast reproduced quickly in ripened 
tea samples, indicating that most of the macromolecules in the tea ma- 
trix were hydrolyzed or converted into small molecules during the 
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3.3. Comparison of bacterial community and richness during the 
manufacturing process and storage 


Another main microorganism in Pu-erh tea is bacteria. The quality- 
filtered bacterial sequences of all eight Pu-erh tea samples were clus- 
tered into 604 bacterial OTUs at the 97%-similarity cut-off. The Venn 
diagram analysis showed that the eight Pu-erh tea samples shared 24 
bacterial OTUs and the bacteria were increased during fermentation 
processes (Fig. 2A). Moreover, the lowest OTU number, 78, was found in 
RMT sample. The Chao index suggested the highest bacterial richness 
was in the RT180d sample (Fig. 2B). The Shannon indices indicated that 
the RT180d sample displayed the highest bacterial diversity (Fig. 2C). 
The bacterial rarefaction curves fully approached the saturation plateau 
(Fig. S1B). The Good’s coverage estimated for the bacterial diversity in 
all Pu-erh samples was higher than 99% suggesting that a major part of 
the diversity in all samples had been captured. 

Contrary to the fungal community, the bacterial richness and di- 
versity increased with extending fermentation times. The predominant 
bacteria genus differed during the fermentation processing and ripened 
(Fig. 2D). A total 604 bacterial OTUs were identified to be assigned into 
21 different phyla, 38 classes, 96 orders, 178 families, 334 genera and 
456 species. Firmicutes, Cyanobacteria, Actinobacteria and other uncate- 
gorized bacteria were the predominant phyla in all samples. Cyanobac- 
teria was the dominant phylum in the raw Pu-erh samples. Proteobacteria 
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and Cyanobacteria were co-dominant in the PF2nd and PF3rd samples, 
while Firmicutes became dominant in the PF4th and PF5th samples 
(>90%). Firmicutes and Actinobacteria were most abundant in the 
ripened tea samples. 

At the genus level, norank_o_Chloroplast was the absolute dominant 
bacteria in the raw Pu-erh tea sample. The norank_o_Chloroplast was 
sharply decreased after entering the manufacturing process, and instead, 
Pantoea abundant increased to 56.68% in PF3rd period. Most bacterial 
richness was decreased in the PF4th and PF5th stages (Fig. 2E), but in 
contrast Bacillus dramatically reproduced to 82.97% in the later stage of 
the fermentation process. The dry and suitable temperature promoted 
bacterial diversity in the ripened period. Staphylococcus and Kocuria had 
synergistic effects on growth. The nutritional resource from the digested 
tea matrix contributed to bacterial community differences between the 
raw and manufacture samples. 

Differences in the microbial community structures among Pu-erh tea 
samples were assessed using Bray-Curtis distance. Based on ANOSIM/ 
Adonis dissimilarity analyses (Figs. S2A and B), the sample distance 
distribution within groups was smaller than between groups. The results 
indicated that all the experimental data had statistical significance. 
Moreover, the RMT sample exhibited the highest fungal diversity but 
lower richness. It is noteworthy that the lowest fungal diversity 
appeared in the fermentation stage. Compared to fungal diversity, bac- 
terial diversity was much higher in the whole manufacturing process. 


Fig. 2. Comparisons of bacterial composition 
and richness in different samples of Pu-erh tea. 
(A) Venn diagram analysis of bacterial OTU 
level in samples during manufacturing process 
of Pu-erh tea. (B) Bacterial taxonomic compo- 
sitions showing the bacterial successions at 
genus level. The bacterial communities in eight 
period samples are comparably diverse. Chao 
index (C) and Shannon index (D). Heat map 
analysis of bacterial (E) communities in Pu-erh 
tea samples at the species level. (note: The 
taxonomic abundance <1% were classified in to 
“others”; “*” P < 0.05 ; “**” P < 0.01; “***” P 
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The result demonstrated that both the diversity and richness of bacteria 
were significant increased in the end manufacturing process. 


3.4. qPCR detection of fungi and bacteria during processing and storage 
of Pu-erh tea 


To further confirm the comprehensive evolution of microorganism 
richness during the processing and storage of Pu-erh tea, we utilized 
qPCR to directly quantify copy numbers of microbial genome for 
different samples. First, we constructed standard plasmids for genera- 
tion of fungal or bacterial standard curves. The p416-TEF-ITS(188) and 
p416-TEF-16S rDNA(659) were constructed for generation of fungal and 
bacterial standard curves, respectively. The fungal standard curve 
(Fig. 3A) and bacterial standard curve (Fig. 3B) were generated by cycle 
threshold (Ct) value. Using genomic DNA extracted from the different 
Pu-erh tea samples as templates, quantitative PCR was performed with 
the same primer sets as those of the standard curve qPCR. The results 
exhibited that fungi amounts increased in the middle of the pile 
fermentation period and then decreased during storage (Fig. 3C), but 
bacteria amounts increased significantly in the initial stage, subse- 
quently decreased and finally increased (Fig. 3D). The 10°?° copies/pL 
of fungi in RMT was significantly higher than the 10°°° copies/pL of 
bacteria. In the PF4th period, the fungal numbers reached the highest 
point at 10°* copies/jiL (Fig. 3C). The numbers of bacteria displayed 
two peaks, 10°7° copies/jiL in the PF2nd and 10?*® copies/uL in the 
RT360d (Fig. 3D). These results were consistent with the overall trends 
of fungi and bacteria revealed by high-throughput sequencing. Our re- 
sults indicated that fungi and bacteria were differently distributed in 
processing stages and worked together to complete the fermentation 
process. 


3.5. Predictive metabolic pathways based on microbial compositions in 
Pu-erh tea samples 


The biochemical composition of Pu-erh tea is mainly determined by 
the types and quantities of microorganisms. After understanding the 
dynamic evolution of Pu-erh tea microorganisms, we predicted and 
analyzed metabolic pathways of microorganisms. The typing analysis at 
OTU level indicated that the samples could be classified into three 
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cluster according to PC1 and PC2 axes (Fig. 4A). Fungal flora classified 
into nine types and community structures were shifted three times 
during manufacturing process (Fig. 4A). 

From the PCA score plot, the samples were divided into four suc- 
cessive processes and eight bacterial types were identified by total OTU 
numbers (Fig. 4B). The result revealed that the bacterial community 
structures were gradually shifted, but enormous changes occurred at 
later fermentation stages and possessed distinct separation with previ- 
ous tea samples. In addition, we found the relativity of bacterial com- 
munity structures between PF5th tea sample and ripened tea sample. 

Fungal communities participate in hundreds of metabolic activities 
throughout the manufacturing process which closely relate to the for- 
mation of special flavors of Pu-erh tea. Through PICRUSt2 (http://hutt 
enhower.sph. harvard.edu/galaxy) function prediction (Langille et al., 
2013), the metabolic activity that the fungal community participates in 
can be effectively obtained, and the longitudinal column of each heat 
map represents the number of various metabolic pathways in the KEGG 
database. Each number corresponds to the detailed description of the 
metabolic pathway in Table S3. The correlation between metabolic ac- 
tivities and samples was manifested by the color depth (Fig. 4C). We 
found that more than 90% of metabolic activities occur in the PF5th 
period and proved the pile-fermentation process played an important 
role in the quality formation of Pu-erh tea. In contrast, the metabolic 
activity of fungi exhibited inactivity in the RMT and ripened period due 
to the storage condition being not appropriate for growth. Overall, the 
metabolic activities of fungi included aerobic respiration I and II, fatty 
acid, beta-oxidation, palmitate biosynthesis I, guanosine nucleotides 
degradation II, t-valine biosynthesis, chitin degradation to ethanol and 
p-galactose degradation V. 

On the other hand, the different KEGG pathways predicted for each 
bacterial community are shown in Fig. 4D. Within the metabolism 
category, the relative abundance is shown of pathways associated with 
energy production and conversion, RNA processing and modification, 
cell cycle control, cell division, chromosome partitioning, lipid transport 
and metabolism, amino transport and metabolism, nucleotide transport 
and metabolism, carbohydrate transport and metabolism, coenzyme 
transport and metabolism. It is worth noting that fungal metabolism is 
mainly involved in the hydrolysis of macromolecular compounds, while 
bacterial metabolism carried out decomposition, transformation and 
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Fig. 3. Quantitative detections of fungi and bacteria performed by qPCR during processing and storage of Pu-erh tea. (A) Standard curve of fungal plasmid; (B) 
Standard curve of bacterial plasmid; Quantitative detections of fungi (C) and bacteria (D) during processing and storage of Pu-erh tea. 
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Fig. 4. Predicted metabolic pathways based on microbial compositions in Pu-erh tea samples. Classification of fungal (A) and bacterial (B) flora based on OTU level. 
Heat map of fungi involved in metabolic pathway (C). Statistical chart of bacteria participating in metabolic activities (D). 
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transfer of small molecular compounds in the Pu-erh tea manufacturing 
process. 


3.6. Metabolomics analysis of Pu-erh tea samples during processing and 
storage 


The color and flavor of Pu-erh tea from different periods displayed 
significant difference, which implied biochemical compositions were 
different among these samples (Fig. 5A). To investigate the metab- 
olomics changes during Pu-erh tea processing, UPLC-MS based metab- 
olomics analysis coupled with multivariate analysis was applied and 
identified critical metabolites responsible for metabolomics variation 
caused by Pu-erh tea processing. The typical total ion current chro- 
matogram for each tea sample is shown in Fig. 5B. Metabolites in tea 
samples were identified according to authentic standards and a tea 
metabolomics database. A total of 209 metabolites were identified in 
eight samples (Table S4). A Venn plot constructed using these data 
suggested that 108 of 209 metabolites were detected in all tea samples 
(Fig. 5C). Among the remaining 101 metabolites, 27 metabolites were 
specific to the RMT sample. These results indicated that Pu-erh tea 
processing caused qualitative changes in metabolic profiles. 

In order to obtain a general overview of the variance of metabolites 
in this study, the 209 metabolites were subjected to multivariate anal- 
ysis using SIMCA-P14.1 ( Umetrics AB ) statistical software. The hier- 
archical cluster analysis (HCA) gave a clear differentiation between raw 
material tea groups and fermented tea groups (Fig. 5D). The discrimi- 
nability revealed by the HCA plot also verified the predictive accuracy of 
the principle component analysis (PCA) and orthogonal partial least 
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squares discriminant analysis (OPLS-DA) models. In brief, the HCA 
grouped the eight samples into two clusters. The raw material tea 
samples were visually different from the other samples, and the PCA and 
OPLS-DA (Fig. 5E and F) analyses also showed a similar clustering. The 
OPLS-DA comparison between the raw material and the other 7 groups 
suggested obvious metabolic differences between the classes in each 
component. The metabolite with a variable importance in projection 
(VIP) value was greater than 1 illustrates a significant contribution to 
the distinguishing of groups within PLS-DA models (Szeto, Reinke, 
Sykes, & Lemire, 2010). The VIP plots revealed that some of the iden- 
tified metabolites contributed to the group separation (Fig. S3). On the 
basis of OPLS-DA results with good pairwise discriminations and the 
cutoff VIP value (VIP>1), a total of 52 paired retention time-mass to 
charge ratio (RT-M/Z) metabolites were selected (Table S5). 

The major metabolic perturbations caused these discriminations 
consist of 52 components. Among these identified metabolites, 12 
catechin, 22 flavonoids and flavonosides, 14 phenolic compounds and 4 
uncategorized compounds. Furthermore, 50 compounds were listed in 
the heat map to provide a more immediate view of the critical metab- 
olites (Fig. 5H). The content variations in the critical metabolites among 
the samples were represented by the depth of the color. The contents of 
key compounds were broadly elevated in the raw material tea samples, 
with exceptions of catechin, myricetin, phloroglucinol, pyrogallol and 
salicylic acid. Compared with raw material tea, most of the key com- 
pounds in Pu-erh tea decreased significantly during the fermentation 
and ripened period. The content of catechins and gallocatechin 
increased slightly during post-fermentation but markedly decreased 
after the middle-fermentation and ripened period. 
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Fig. 5. Metabolomics analysis of Pu-erh tea samples during processing and storage. (A) The pictures of dried tea materials and extracts from different periods. (B) 
Total ion chromatograms of Pu-erh tea samples in different periods. (C) Venn plot of the identified metabolic compounds in different tea samples. (D) Multivariate 
analysis of eight samples, dendrogram plot of tea samples. (E) PCA score plot. (F) OPLS-DA score plot. (G) Permutation plot of OPLS-DA. (H) Heatmap analysis of 
critical metabolites in different samples. Each column represents a tea sample, and each row represents a critical metabolite. A color-coded scale grading from green 
to red corresponds to the content of critical metabolite shifting from low to high. (For interpretation of the references to color in this figure legend, the reader is 


referred to the Web version of this article.) 
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Fig. 5. (continued). 


3.7. Biochemical measurements of critical functional components of Pu- 
erh tea 


UPLC-MS approach mainly focuses on the identification of small 
molecular endogenous metabolites in organisms. However, the nutri- 
tional and healthy ingredients in Pu-erh tea are mostly macromolecules. 
To further confirm critical functional components of Pu-erh tea, 
biochemical methods were used to measure the contents of these 
components. 


3.7.1. Polyphenol 

Phenolic compounds are important aromatic metabolites contrib- 
uting to the flavor and anti-obesity activity of Pu-erh tea (Oi, Hou, 
Fujita, & Yazawa, 2012). Polyphenols (218.86 mg/g) in raw materials 
sharply decreased to 73.76 mg/g in RT360d samples (Fig. 6A). These 
results might be caused by oxidation of tea microorganisms. Tea 


polyphenols are decomposed into chatein, chalignin and phenolic acid. 
On the other hand, they are oxidized into polyphenols complex. Cate- 
chins account for 60-80% of tea polyphenols and are the main compo- 
nents contributing to the antioxidant activity of tea (Chen, Liu, & Chang, 
2010). During dark tea processing, polymeric catechins and gallyolated 
catechins are gradually hydrolyzed into non-gallyolated catechins and 
gallic acid, which are further degraded into phenolic acids (Ge et al., 
2019; Li, Dai, Xie, & Tan, 2018). Notably, phloroglucinol and pyrogallol 
are generated during later fermentation and may have originated from 
the degradation of catechins, resulting in significantly changed poly- 
phenol during the manufacturing process. 


3.7.2. Amino acid and tea protein 

Amino acids are important ingredients contributing to the taste of tea 
infusion. Additionally, y-Aminobutyric acid, a non-protein amino acid, 
acts as a major inhibitory neurotransmitter in microbial fermentation 
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Fig. 6. Biochemical measurements of functional components of Pu-erh tea. (A) 
Changes of tea protein, total flavonoids, tea polyphenols, tea polysaccharide 
and free amino acid and (B) tea pigment (theaflavins, thearubigins, thea- 
brownins) contents during processing and storage of Pu-erh tea; the data in 
three replicas, (a-f ) expressed P < 0.05. 


(Zhao et al., 2011). The time-varying change of 21 amino acids was 
shown in (Table S6). Total amino acids were decreased from 68.94 mg/g 
in the RTM sample to 16.59 mg/g in the final ripened samples (Fig. 6A). 
Previous study found that threonine and aspartic acid were positively 
correlated with the umami taste of the tea infusion (Yang et al., 2018). 
Thus, the amino acid decrease may account for the reduction of umami 
and fresh tastes in Pu-erh tea infusion. The amino acid reduction could 
be due to decarboxylation, deaminization, other reactions forming new 
substances, or consumption as nitrogen sources by microorganisms. 
Contrary to the decrease of amino acids, the contents of tea protein 
increased significantly during the storage period. 


3.7.3. Flavonoids and flavone glycosides 

The total contents of flavonoids were markedly reduced during Pu- 
erh tea processing (Fig. 6A). This may be caused by the hydrolysis of 
flavonoids and flavone glycosides into catechins and gallic acid (Qin, Li, 
Tu, Ma, & Zhang, 2012). The flavonoids and their derivatives may be 
converted into microbial metabolites during microbial fermentation 
(Zhang et al., 2011). A previous study found that flavan-3-ols and pro- 
cyanidin B3 contributed to the bitter and astringent tastes of tea infusion 
(Yang et al., 2018). The reduction of flavonoids and flavone glycosides 
may improve the bitter and astringent tastes of Pu-erh tea. 


3.7.4. Water-soluble polysaccharide 

Water-soluble tea polysaccharide is one of the main compounds 
contributing to the tea flavors that can relieve the bitter taste from 
polyphenols, theobromine, chalignin and other ingredients. The con- 
tents of tea polysaccharides were gradually increased in the post- 
fermentation stage (Fig. 6A). The reason was that tea polysaccharides 
in the RMT are mostly combined with protein and lipid to form water 
insoluble complex sugar which is digested by hydrolytic enzymes from 
Aspergillus niger, or due to the large reproduction of fungi in early 
fermentation, as their cell walls contain many polysaccharide com- 
pounds. As the fermentation proceeded, the contents of the water- 
soluble polysaccharides decreased to 32.74 mg/g in the RT360d sam- 
ple. We noticed that the content of sucrose was higher in raw tea owing 
to soluble sugar being utilized by the microorganism as a major source of 
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carbon during the process of fermentation. 


3.7.5. Other notable changes 

Tea pigments are complex and group compounds, which are difficult 
to separate into single compounds. So the analysis of main tea pigments 
including theabrownin (TB), theaflavin (TF) and thearubigin (TR) was 
performed with a system approach (Xie et al., 2009). The contents of TF 
and TR were significantly decreased, both of which contribute to the 
astringency and bitterness of tea (Fig. 6B). On the contrary, the increase 
of tea pigment TB gives a characteristic brown color of Pu-erh tea and 
plays an important role in the quality of the final product. 


3.8. Correlation between microorganisms and metabolites during 
manufacturing process and storage of Pu-erh tea 


The flavor and quality of Pu-erh tea depends on metabolites deter- 
mined by compositions of microorganisms. In this study, the most 
important target is to fully explore the correlation between microor- 
ganisms and metabolites during manufacturing process and storage of 
Pu-erh tea. The critical metabolites, including tea polysaccharide (TPS), 
theobromine, caffeine, protein, TBs, catechins, flavonoids, amino acids 
(Aa), polyphenols, theanine and TRs were analyzed in the top50 fungi 
and top50 bacterial (Fig. 7A). The correlation between Cladosporium and 
GCG metabolism reached 90% (Fig. 7B). Aspergillus niger plays a critical 
role associated with the metabolism of almost all macromolecular 
compounds, such as polyphenols, flavonoids, pigments, theanine and tea 
polysaccharides in the post-fermentation period (Fig. 7B). This suggests 
that Aspergillus niger played a critical role in the formation of Pu-erh tea 
quality in the early stage of the fermentation process. Rasamsonia_e- 
mersonii is considered as the principal chalignin producing microor- 
ganism during the later fermentation process. There was a positive 
correlation between the abundance of Blastobotrys adeninivorans and 
Thermomyces_lanuginosus to the contents of alkaloid (caffeine, theobro- 
mine) and the large accumulation of tea protein. In our study, the effect 
of fungi on 19 metabolites was clarified and each metabolite corre- 
sponded to one or more fungus. 

Bacteria also affect the metabolism of endogenous compounds in tea. 
The metabolic activity (80%) is associated with dominant bacteria in 
RMT and early fermentation stage (Fig. 7C). Cyanobacteria and Pantoea 
participated in the metabolism of macromolecular compounds such as 
tea polysaccharides, flavonoids, theanine, polyphenols etc. The domi- 
nant Bacillu was mainly responsible for gallic acid and theobromine 
production in the late fermentation stage of the Pu-erh tea 
manufacturing processes. Staphylococcus and Kocuria are involved in the 
accumulations of tea protein, chalignin and caffeine in the ripened 
period (Fig. 7D). Moreover, the effect of bacteria on 19 metabolites is 
not strong as that of fungi, and they may play an auxiliary role in sub- 
stance transformation. For example, Bacillus and Staphylococcus pro- 
duced some oxidases which similarly accelerated the transformation 
reaction and reduced the fermentation cycle. 


3.9. Correlation between environmental factors and microorganisms 


Environmental conditions regulate compositions and abundance of 
microorganisms, thereby affecting the flavor and quality of Pu-erh tea. It 
is critical to explore the correlation between environmental factors such 
as temperature (Fig. 8A), humidity (Fig. 8B), pH (Fig. 8C), and moisture 
(Fig. 8D) and the dominant fungi and bacteria in each sample. The result 
shown that the growth of Cladosporium in RMT was basically unaffected 
by these environmental factors. During the processing and storage of Pu- 
erh tea, both dominant of fungus Aspergillus niger and Rasamsonia_e- 
mersonii were regulated by temperature and humidity (Fig. 8E). The 
optimal condition for reproduction of these two species was 50-55 °C 
and humidity of 30%-35%. The dominant Blastobotrys_adeninivorans 
was regulated by pH, which displayed a low correlation with tempera- 
ture and humidity. In addition, nearly 90% of the fungi were regulated 
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Fig. 7. Correlation between microorganisms and metabolites during process and storage of Pu-erh tea. (A) PCA analysis and (B) heat map analysis of the correlation 
between 19 critical metabolites and fungi at the species level. (C) PCA analysis and (D) heat map analysis of the correlation between 19 critical metabolites and 


bacterial on genus level. 


by the environment, but less than 10% of dominant fungi manifested a 
strong correlation with temperature and humidity (Fig. 8E). This also 
explained that it is necessary not only to maintain a suitable temperature 
and humid condition, but also to monitor the pH during the fermenta- 
tion process. pH 5 was the appropriate fermentation condition. 

The dominant bacteria in RMT was Cyanobacteria, which was un- 
regulated by these three environmental factors. Humidity showed the 
greatest impact on Pantoea growth and reproduction in the middle stage 
of fermentation. The dominant bacterium in late fermentation period 
was Bacillu and its optimal growth conditions have certain correlations 
with both temperature and pH. The growth of Staphylococcus and 
Kocuria in the ripened period were largely associated with pH. In sum- 
mary, pH, temperature, and humidity exhibited a great effect on the 
growth and reproduction of bacteria during the processing and repined 
periods (Fig. 8F). Therefore, to increase the quality and enhance the 
growth of beneficial bacteria (such as Lactobacillus and Bacillus) during 
fermentation, a suggested temperature of 40-50 °C and pH 5 could be 
applied in manufacturing process (Fig. 8F). 


3.10. Correlation between microorganisms and bioactivity of Pu-erh tea 


To explore the antioxidant ability of Pu-erh tea samples in different 
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periods, the DPPH (Fig. S4A) and ABTS (Fig. S4B) free radical scav- 
enging ability of in different samples were conducted. In addition, we 
analyzed the correlation of DPPH and ABTS radical clearance with the 
top 6 most abundant fungi and bacteria. The results showed (Fig. S5) 
that the highest correlation with DPPH radical clearance was Blastobo- 
trys_adeninivorans and the correlation coefficient was 0.8675 (Fig. S5D). 
Next were Thermomyces_ lanuginosus and Aspergillus fumigatus, with 
correlation coefficients of 0.5429 (Fig. S5E) and 0.5369 (Fig. S5C), 
respectively. The lowest correlation was Rasamsonia_emersonii 
(Fig. SSA). Similarly, the highest association with ABTS radical clear- 
ance was also Blastobotrys_adeninivorans with correlation coefficient of 
0.7164 (Fig. S5D), followed by Thermomyces lanuginosus with correla- 
tion coefficient of 0.6116 (Fig. S5E), and the remaining four fungi 
showed less correlation. Therefore, Blastobotrys.adeninivorans and 
Thermomyces_lanuginosus are the key strains to enhance the antioxidant 
resistance of Pu-erh tea. The abundance proportion of these two strains 
could be appropriately improved in the fermentation process. The cor- 
relation analysis between antioxidant capacity and relative abundance 
of top 6 bacteria indicated (Fig. S6) that the removal of these two free 
radicals was not in a high correlation with bacteria. The possible reason 
is that bacteria rarely undertake the transformation function of macro- 
molecular compounds in the fermentation system. Instead, bacteria are 
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Fig. 8. Correlation between environmental factors and microorganisms. Change of the temperature (A), relative humidity (B), pH (C) and tea moisture content (D) 
during the manufacturing process of Pu-erh tea. Heat map analysis of the correlation between environmental factors and abundance of top50 Fungi (E) and bac- 


teria (F). 


responsible for decomposition of small molecules and production of 
enzymes which catalyze the transformation reaction leading to reducing 
the fermentation cycle. Therefore, the correlation between bacteria and 
antioxidants is very low. 

We analyzed the correlation between DPPH and ABTS radical 
clearance with the six main active ingredients. It reported that the main 
active substances associated with antioxidant activity include poly- 
phenols, theabrownins, catechin, flavonoids, polysaccharides and free 
amino acids. The results showed (Fig. S7) that the highest correlation 
with DPPH radical clearance was free amino acids and the correlation 
coefficient was 0.7465 (Fig. S7F) followed by polysaccharides with 
correlation coefficient of 0.4932 (Fig. S7E). The correlation of poly- 
phenols, flavonoids and catechins on DPPH radical clearance were 
almost similar. Chalignin had the lowest correlation coefficient of 0.3 
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(Fig. S7C). The polyphenols displayed the highest correlation with ABTS 
radical clearance with correlation coefficient of 0.8623 (Fig. S7A) fol- 
lowed by free amino acids with correlation coefficient 0.8188 (Fig. S7F). 
Chalignin and flavonoids had similar ABTS radical clearance. The lowest 
correlation coefficient of polysaccharide was only 0.5856 (Fig. S7E). In 
summary, the ingredients most associated with antioxidant properties 
were tea polyphenols and free amino acids, but both are decreased 
during fermentation stage. We suggest a new avenue to improve the 
quality of Pu-erh tea by increasing the content of polyphenols and free 
amino acids during the manufacturing process. 


4. Conclusion 


Unravelling the microbial communities responsible for producing 
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beneficial bioactive compounds is critical in Pu-erh tea manufacture. 
Considering the lack of a systematic comparative study during the Pu- 
erh tea manufacturing process, we carried out a comprehensive multi- 
variate analysis focusing on microbial community structure shift, cor- 
relation between microbiomes and metabolites, and antioxidant 
changes during Pu-erh tea processing. The raw material tea sample 
displayed the highest fungal diversity, while lower richness and the 
lowest fungal diversity appeared in fermented tea. Compared to fungi, 
the diversity of bacteria was much higher across the whole 
manufacturing process. Both the diversity and richness of bacteria were 
significantly increased at the end of the manufacturing process. The 
conversion of metabolites in Pu-erh tea was mainly catalyzed by mi- 
crobial enzymes secreted from microorganisms during the 
manufacturing process. A total of 50 critical metabolites were respon- 
sible for metabolic changes caused by Pu-erh tea processing. The con- 
tents of polyphenols, polysaccharides and flavonoids significantly 
decreased while theabrownin, some novel phenolic acids and catechin 
derivatives formed during processing. The environmental condition 
displayed a great effect on the growth and reproduction of bacteria, but 
less than 10% of dominant fungi exhibited strong correlation with 
temperature and humidity. In summary, our results reveal the microbial 
composition is the critical factor in changing the metabolic profile and 
provide a theoretical basis for improving the quality and safety of Pu-erh 
tea. 
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